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This paper provides a theoretical investigation on the thermodynamic analysis of a Stirling engine with 
linear and sinusoidal variations of the volume. The regenerator in a Stirling engine is an internal heat 
exchanger allowing to reach high efficiency. We used an isothermal model to analyse the net work and 
the heat stored in the regenerator during a complete cycle. We show that the engine efficiency with 
perfect regeneration doesn't depend on the regenerator dead volume but this dead volume strongly 
amplifies the imperfect regeneration effect. An analytical expression to estimate the improvement due to 
the regenerator has been proposed including the combined effects of dead volume and imperfect 
regeneration. This could be used at the very preliminary stage of the engine design process. 

© 2010 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine was invented in 1816 by Robert Stirling [1], 
The first era of the Stirling engine was ended by the rapid devel¬ 
opment of internal-combustion engine and the electric motor. The 
second era began around 1937 [2j. A wide range of investigations 
has been made on the development of Stirling engine [3] but 
further developments are still needed for wider practice applica¬ 
tions [4—6], Its application has several advantages like multi-fuel 
capability, low fuel consumption and low noise [3], The variety of 
fuel sources (gases or liquids) and abundant heat sources (solar 
radiation [7]) that can be used are of particular interest. The current 
applications under development are heating and cooling, combined 
heat and power, solar power generation, Stirling cryocoolers, heat 
pump, marine engines, nuclear power generation, automotive 
engines, electric vehicles, aircraft engines and low temperature 
difference engines. Commercial Dish-Stirling system of 25 kW [8] 
takes the energy from the sun and converts it into electrical 
energy providing a clean, efficient and reliable source of energy. 
Micro systems are of primary importance in the latest develop¬ 
ment. Demonstration models have been realized like commercial 
Low Temperature Difference Stirling Engines running on any warm 
surface (i.e„ a cup of hot water) [9], Reduction in the dimensions 
can lead to new designs. The moving component could be replaced 
by vibrating membranes as in the miniaturized membrane engine 
[ 10 ]. 

Despite the apparent conceptual simplicity of the Stirling 
engine, the mathematical analysis is complicated. The first analysis 
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of a Stirling engine was published by Schmidt in 1871 [11] and has 
been reported in a more comprehensive form in the book of Urieli 
and Berchowitz [12]. For all new designs, a simple thermodynamic 
analysis is performed even in the case of recent development to 
determine the characteristics [13]. As efforts to reduce some loss 
mechanisms could tend to aggravate others, complementary 
approaches are needed. A simple model could be used at the first 
stage followed by a sophisticated approach including computation 
could be necessary. It is noted that a practical Stirling cycle always 
departs significantly from the ideal cycle mainly due to the 
performance of the heat exchanger and the regenerator. Some 
detailed analysis with practical values and examples of the engine 
are given by Organ [14]. The commercially available codes, as CFD- 
ACE and Fluent, provide a macroscopic porous-media model where 
permeability, inertial coefficient, and porosity are required as the 
input parameters. Pressure drop and temperature profile are 
obtained with accuracy [15] and could be used to improve the 
existing system [16] Anyway, a simple mathematical model calcu¬ 
lation is of primary importance to understand in depth the 
parameters that need to be optimized. 

Several forms of cylinder coupling in Stirling engine can be used, 
but the more common is the opposed-piston configuration known 
as alpha configuration with sinusoidal volume variation of working 
space [3[. 

Here we consider a model system of practical importance con¬ 
sisting of a notional discontinuous piston motion and a simple 
harmonic piston motion. This work could be used in the under¬ 
standing of an imperfect regenerator and the dead volume coupled 
effects in the primary stage of an engine design. It is a comple¬ 
mentary work of Kongtragool and Wongwises [17] who have 
associated the heat of the hot isothermal expansion, a part of the 
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Nomenclature 

8 

cycle angle varying between 0 to 2n 



T 

temperature 

yi 

surface section of the regenerator 

T h 

heater temperature 

Cv 

specific heat at constant volume 

Tk 

cooler temperature 

Cv, m 

molar specific heat at constant volume 

T r 

regenerator temperature 

e 

thermal efficiency 

V 

volume 

T 

ratio of the specific heat at constant pressure to 

v h 

heater volume 


constant volume 

V k 

cooler volume 

L 

length of the regenerator 

V m 

minimum volume considering V h and V k only 

n 

mole number 

Vm 

maximum volume considering V h and V k only 

n h 

mole number in the heater 

V r 

regenerator volume 

n k 

mole number in the cooler 

W 

work 

n r 

mole number in the regenerator 

Wh 

work for the heater piston 

P 

absolute pressure 

w k 

work for the cooler piston 

On 

heat furnished by the heater 

w r 

engine net work 

Qk 

heat furnished by the cooler 

X 

fraction of the heat furnished by the heater and 

Or 1-.3 

heat released by the regenerator 


released by the regenerator 

R 

gas constant 

Z 

linear abscissa 

U 

internal energy 




cycle, to the heat supplied externally in order to explore the 
combination of dead volume and imperfect regeneration effects. 
However, the mole variation in the regenerator during isothermal 
transformations has not been properly taken into account. In this 
article, we have carefully analysed the combination of dead volume 
and imperfect regeneration particularly the effect to the heat stored 
in the regenerator and the mole number in each of the three 
volumes (heater, cooler and regenerator) that enables us to 
understand the variation of efficiency in these cases. With this 
simple approach, we compare the heat stored and the heat rejected 
by the regenerator with the net work. This allows us to determine 
an analytical expression that takes into account the dead volume 
and the imperfect regeneration. 

2. Stirling cycle analysis 

The Stirling cycle based on two opposing pistons in motion is 
shown in Fig. 1. It consists of two isothermal and two isochoric 
processes (Fig. 2). During the isothermal compression (1—2), the 
cold piston moves towards the regenerator. In the following step 
(2—3), both pistons move similarly so that the volume remains 
constant and the working fluid is transferred through the regen¬ 
erator. This step is followed by the isothermal expansion process 
(3—4) with only the hot piston moving. In the step (4-1) both 
pistons move keeping the volume constant and transferring the 


Regenerator 




working fluid through the regenerator. This is an ideal cycle 
assuming the perfect heat transfer and does not take into account 
heat losses during the regeneration. First we consider an ideal case 
of perfect gas with infinite conductivity for all materials. Therefore, 
we have: 

dU = C v dT = nC v m dT, <5 W = -pdV and 5Q 
= dU - 6W 

2.3. Regenerator temperature 

Within the isothermal hypothesis, the temperature should vary 
continuously between the cold ( T k ) and the hot space (Tti). Irre¬ 
versible processes take place if a temperature difference exists 
between the thermostat and the gas. Here, the gas temperature is 
always at the thermostat temperature, leading to no creation of 
entropy. Experimentally, for homogeneous regeneration, a linear 
variation of the temperature versus the position has been observed 
[12]. Considering that the temperature varies linearly from T k at the 
abscissa 0 to Tj, at the abscissa L, the mole number n r (z) in the 
volume Adz at the abscissa z is: 

n (z) = pAdZ 
r{> R(T h +f(T k - T h )) 

In the regenerator volume V — AL , n r is 



Fig. 2. Watt Stirling cycle with notional discontinuous piston motion. 
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"'-Jw-rn wi,h 

Consequently, we can use the regenerator temperature for all 
mechanical calculations, i.e., to determine pressure. If the hot and 
the cold space are at the same temperature, the divergence in the 
expression can be avoided by adding an incremental value to one of 
the temperatures. For heat exchange, with no dead volume, all 
molecules go through the regenerator and receive an amount of 
heat dependent only on T k - T k and the regenerator temperature 
has no effect. Nevertheless, the direction of the heat flow is dictated 
by the second law of thermodynamics. 

3. Notional discontinuous piston motion 

3.3. Without regenerator dead volume and with perfect 
regeneration 

In the simplest case without dead volume and with perfect 
regeneration, all quantities can be easily calculated (Table 1 ). 
During the isochore transformation, the pressure is given by 

nR 

P ~ 

T h T k 

So, we can calculate the work W k = f pdV k and W k = J pdV k 
which are opposite in sign to each other (see Table 2) 

The Carnot efficiency is reached according to: 

e = ~ = 1 - ^ with W t = W h + W k 
til l h 

= -"<«>- 7 .)'n(g) 

where W t is the net work and e the efficiency. 

3.2. With regenerator dead volume and with perfect regeneration 
In this case, the mole number is given by: 


n = n k + n r + n h , 


Table 2 

Work for cold and hot space without dead volume. The work for a transition is 
simply W = W h + W k . 


Transition 

Work hot space W/, 

Work cold space W k 

1 -*2 

0 


2—>3 

-nR^Jr 

•r 

nR 7 ^ 

b 

3-*4 


0 

4—►! 

nR^ 

-nR 7 ^ 

b 


To determine the heat transferred by the hot space, the 
isothermal transition 3—>4 cannot be used as the mole number in 
the regenerator is not identical at position 3 and position 4. 
However the cyclic change in the number of moles of the working 
gas has to be zero for each of the cells. Thus, integrating over the 
cycle, we obtain for the working spaces: 

Qk = -W k , Q/, = -W h and Qr = 0 


For the ideal regenerator Q, = 0 because all heat exchange 
between the regenerator matrix and the working gas is internal. A 
first order development of the net work versus the regenerator 
dead volume gives: 


W t 


-nR(T h - T k ) 



-V r 


Jh + T k 

T r 




It is observed that the net work is reduced. 

At the same time, the thermal efficiency equals the Carnot 
efficiency 


W t _ W h 34 + W k i2 _ 1 T k 
Oh ~ W h u + W h23 + W M 4 - T h 

So, we can conclude that the dead volume has no influence on the 
efficiency in the ideal case, that contradicts to the results obtained 
by Kongtragool and Wongwises [17] where a decrease of the effi¬ 
ciency is found with increasing dead volume. The disagreement 
comes from the simplifications used by Kongtragool and Wongw¬ 
ises. Contrary to our approach, the authors [17] consider a constant 
mole number in the regenerator during the isothermal expansion. 


the pressure by 
n 

P = Yl.Yl . V 

RT k RT r RT h 
and the internal energy by 

U = —^ An k T k + n r T r + n h T h ) with T = ^ 

7-1 <-v 

The piston work can be calculated and is presented in Table 3. 


3.3. Without regenerator dead volume and with imperfect 
regeneration 

Taking into account all the conclusions from the previous 
subsection, one can write the added heat as: 

Qr,1^3 = U 3 -U 2 = nC V m (T h - T k ) = n^(T h - T k ) 

Table 3 

Work for cold and hot spaces with dead volume. 


Energy quantities without dead volume. 


Transition 

Work W 


Heat Q 

A U 


1 —>2 

-nRT k ln( 

'V m \ 

Wm) 

">'"0 

0 

2—*3 

2—>3 

0 


tlCv, m(Ti i — T k ) 

nCy, m(Th - T k ) 

3-»4 

3—>4 

-nRT k ln( 

'Vm\ 

\Vm) 

nRr ‘ ln © 

0 


4—>1 

0 


nCy, m {T k - Th) 

m(T k — Th) 

4—>1 


Work hot space W h 


0 



-nR^r 


7J - 7* 



Work cold space W k 


-nRT k In 


nR-r 


+ 
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Fig. 3. Heat exchange in the regenerator. 


If x is the fraction of the external heat transferred to the regenerator 
when the gas is heated, we should add this heat to the total heat 
input to the system 

Qinput = *0r,l->3 + Qh 

leading to an efficiency of 


i \ 


e = 


W t 

Qinput 



1 



3.4. With regenerator dead volume and with imperfect 
regeneration 

In order to investigate this case numerical simulations are 
needed. The energy released and stored by the regenerator is 
calculated using V m = 10 L, Vm = 40 L, R = 8, 31, y = 1.4, n = 10 mol, 
Tfc = 300 K, Tf, = 1000 K. The regenerator volume is indicated in the 
graphs in which the regenerator heat variation is plotted versus the 
angle 8 (see Fig. 3). 

For the 3 spaces, we can calculate the variation of the internal 
energy (dU = nCy, m dT), the elementary work <5W = -pdV, and the 
elementary heat as follows: 

• in the hot and cold space: 


5Q = dU-5W , 


With the example reported in Table 4, we clearly see that the 
heat exchange in the regenerator is nearly the same while the net 
work is strongly reduced. So the association of both these effects 
leads to a decrease in the efficiency and the minimal value is 
reduced due to the regenerator dead volume. One can verify that — 
W t IQj , equals the Carnot efficiency 1 — T^/Tf, and satisfies the second 
law Qfi/Tji + QiJTk = 0 with exchanged entropy and no irreversible 
creation. 

In Fig. 4, it has been shown that the mole number in the 
regenerator versus the cycle angle 6 is never constant during 
a transformation. This observation explains why using the 
expression for work is the good method to calculate the efficiency 
and why we need computation to find the heat exchange between 
the gas and the thermostat for the hot, cold and regenerator spaces. 
In Fig. 5, we have reported the efficiency with and without dead 
volume considering imperfect regeneration. The efficiency value at 
x = 0 is the Carnot efficiency. In the case of dead volume, the effi¬ 
ciency decreases strongly with the imperfect regeneration. It is 
surprising to find that with the dead volume in our example, an 
imperfect regeneration of x = 0.28 gives the same efficiency as with 
no dead volume and x = 1 as reported in Fig. 5. 

To calculate the efficiency, we start from 


e = 


W t 

xQr ,1—>3 ' Qh ,cycle 


with Qh .cycle 


leading to 



e = 




1 



• in the regenerator space: 

5Q,- = dU r + dn k ^ r C Vm (T r - T k ) + dn h ^ r C Vm (T h - T r ) 

For V r = 0, we obtain Qr = nC v , m (Th - T/ ( ) = 145.4 kj which is the 
area corresponding to the 2->3 transformation. 


Table 4 

Net work, heat released by the regenerator (half cycle) and heat exchange with 
cooler and heater. 


Vr (L) 

w, (kj) 

Or.l-3 (kj) 

a,(kj) 

Ok (kj) 

0 

-80.4 

145.4 

115.1 

-34.7 

20 

-20.3 

134.9 

29.0 

-8.7 
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Fig. 4. Variation of the mole number versus the cycle angle 0 with and without 

regenerator volume Fig. 6. Effect of the dead volume to the ratio -Q,- i_ 3 /W t . For the approached formula, 

we use the correction reported in equation (1). 


In all cases -Q,-,i-> 3 /W t increases linearly with dead volume (see 
Fig. 6 for an example with the input data given in the beginning of 
the subsection). Consequently, any increase in the dead volume 
amplifies the effect of the imperfection. Calculation of the linear 
coefficient dependencies is complicated because of the small vari¬ 
ation of Or , _ 3 and the large variation of W. If the dead volume is 
not taken into account, we have 


e 



T A _ 

T hJ \+x(V r 


0 ) 


1 


0r,t-.3(Vr = 0) 

-W t (V r = 0) 



than 20%) can be made using W t (V r ) and considering Q, i ^,3 to be 
independent to dead volume, leading to 


X(Vr) 


hi+Ffc 


1 + V T - 


1 

% 


log 


< x(V r = 0) 


( 1 ) 


This approximation is in good agreement with the numerical values 
in Fig. 6. 


and with dead volume 


4. Sinusoidal piston motion 


e 




1 


l+x(Vr) 


0r,l->3 (Vr) 

-W t (V r ) 


0 



So the regenerator with dead volume should satisfy the 
condition 


A Watt Stirling cycle with sinusoidal piston motion is shown in 
Fig. 7. The minimum volume V m and the maximum volume V M are 
also denoted. The sinusoidal piston motions used with a dephasing 
angle of 7r/2 (this value is very close to the angle giving the 
maximum net work) are [11]: 


0r,l-.3 (Vr) 

-W t (V r = 0) 

This condition is of practical importance: the dead volume should 
be as small as possible. A rough approximation (with an error of less 



x 

Fig. 5. Efficiency with and without dead volume and imperfect regeneration. 


v h = Vsw^l + cos(fl)) 

and 
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Fig. 8. Heat exchange in the regenerator. 


V/i — Vsw-^l + COS^0 +2^)) — ^sw^O — cos (^)) 

with 0 varying between 0 and 2n and V sw being the swept volume. 
In this trivial case, the total volume V equals 

v = V sw ^l+^cos(0 + £)j 

leading to 

V m = V sw f 1 —2 j ^rid = V S w 


4.2. With or without regenerator dead volume and with perfect 
regeneration 

Using the approach of Berchowitch and Urieli [12], we have 



c = 


Vsw /J_ ± 

2 vi n 


_ V r Vsw Vsw 
- 7> + 27^ + Z4 


b = 


and 


Pci 


_ 2nR ( i 1 


The heat supplied externally is opposed to the work done by the hot 
piston according to the first law applied over the cycle 
(AU = 0 = Qh + W/,). The Carnot efficiency is reached as follows: 


2n 

W h = - j pdV h = - J P^dd = ^-Vs W pc isin(0) 

cycle 0 

and 


w k = -^VswPc 1COS(0) 
with 


j8 = arctan 



Table 5 

Net work, heat released by the regenerator and heat exchange with cooler and 
heater. 

V r (L) W t (kj) Qr,i -*3 (kl) Qh (kj) Qk (kj) 

0 -68.0 145.4 97.4 -29.4 

20 -17.6 115.5 25.1 -7.6 



Fig. 9. Effect of the dead volume to the ratio - Qr_ >o/W t . 
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e 


w t 

Oh 


Wft - w k 

-w h 


= 1 - tan (3 = 1- 


Tk 

T h 


4.2. With or without regenerator dead volume and with imperfect 
regeneration 


In the case without regenerator dead volume, the analysis is 
identical to the one described above with discontinuous motions of 
the piston. Even if the regenerator dead volume is taken into 
consideration, the same numerical approach can be used. The heat 
exchange in both cases with and without dead volume are pre¬ 
sented in Fig. 8 . 

Table 5 details the energy exchanged during the cycle. With 
dead volume, the absolute values of the work and the heat 
exchanged with the hot and cold thermostats are strongly reduced 
while the heat exchanged by the regenerator is less affected. 

We have used the previous result (equation (1)) and plotted 
a straight line as shown in Fig. 9 corresponding to: 


W t 

Qr,> 0 


<Yr 


0 ) 



Even if the complexity of the problem does not allow us to find the 
exact expression, this linear approximation is an important result of 
this work. To be efficient, the regenerator should have a dead 
volume that is as small as possible. 


5. Conclusions 


A careful analysis has revealed that, in the regenerator with 
a dead volume, the mole number is not constant during the cycle. 
With perfect regeneration, Carnot efficiency is reached irrespective 
of the dead volume but with dead volume, the stored energy in the 
regenerator is high in comparison with the net work. Without 
perfect regeneration, the efficiency reduction due to the imper¬ 
fection is strongly amplified by the dead volume. Using a simple 
model with standard thermodynamic analysis, we have shown that 


a simple expression can be found for characterization of both the 
imperfection (x) and the dead volume (V r ) in the regenerator given 
by the formula: 


x(V r ) 


1 + V r 


_ 1 _ 

T - w m v M 


log 


< x(V r 


0 ) 


This formula is a rough approximation that is useful for discon¬ 
tinuous and sinusoidal motions of the pistons. 
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